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C/} ' Abstract 



We present, using spectral analysis, a possible way to prove the Rie- 
mann's hypothesis (RH) that the only zeroes of the Riemann zeta-function 
are of the form s = 1/2 + i\ n . A supersymmetric quantum mechanical 
model is proposed as an alternative way to prove the Riemann's conjec- 
iq ■ ture, inspired in the Hilbert-Polya proposal; it uses an inverse eigenvalue 

\Q ' approach associated with a system of p-adic harmonic oscillators. An in- 

£SJ , terpretation of the Riemann's fundamental relation Z(s) = Z{\ — s) as a 

[ — ■ duality relation, from one fractal string L to another dual fractal string 

<^^> L' is proposed. 

1 Introduction 

O ' Riemann's outstanding hypothesis (RH) that the non-trivial complex zeroes of 

D , the zeta-function £(s) must be of the form s = 1/2 ± i\ ni remains one of the 

open problems in pure mathematics. The zeta-function has a relation with the 
number of prime numbers less than a given quantity and the zeroes of zeta are 
deeply connected with the distribution of primes |12| . The spectral properties of 
the zeroes are associated with the random statistical fluctuations of the energy 
levels (quantum chaos) of a classical chaotic system J5J. Montgomery Q has 
shown that the two-level correlation function of the distribution of the zeroes 
is the same expression obtained by Dyson using random matrices techniques 
corresponding to a Gaussian unitary ensemble. See also |ll[, [ [18| |24| ] and g. 
An extensive compilation on zeta related papers can be found at |26|. 

One can consider a p-adic stochastic process having an underlying hidden 
Parisi-Sourlas supersymmetry, as the effective motion of a particle in a potential 
which can be expanded in terms of an infinite collection of p-adic harmonic 
oscillators with fundamental (Wick-rotated imaginary) frequencies tu p — i In p 
(p is a prime) and whose harmonics are ujp. n = iln p n (See H). This p-adic 
harmonic oscillator potential allowed to determine a one-to-one correspondence 
between the amplitudes of oscillations a n (and phases) with the imaginary parts 
of the zeroes of zeta X n , after solving a inverse eigenvalue problem. 



Pitkanen m& proposed an strategy for proving the Riemann hypothesis in- 
spired by orthogonality relations between eigenfunctions of a non-Hermitian 
operator that describes superconformal transformations. In his approach the 
states orthogonal to a vacuum state correspond to the zeroes of Riemann zeta. 
However, a proof was not given. 

The contents of this work are the following. In section we consider the 
SUSY-QM model approach to the Hilbert-Polya proposal. In section plwe define 
a different operator than Pitkanen's, expressed in terms of the Jacobi theta 
series, and we find the orthogonality relations among its eigenfunctions, and 
finally we prove some theorems and present a novel approach to prove the RH. 
In section [| a discussion of the fractal string construction, in relation to the 
Riemann zeta-function, given by Lapidus and Frankenhuysen and a possible 
generalization is presented. Finally, some concluding remarks concerning the 
multifractal distribution of prime numbers found by Wolf and the distribution 
of lengths of fractal strings are provided. 

2 A supersymmetric potential 

The Hilbert-Polya idea to prove the RH is the following ||. If the zeroes of 
the Riemann zeta function are s n = 1/2 + i\ n , then it must exist a Hermitic 
operator T such that the s n are complex eigenvalues of the operator 1/2 + iT, 
in other words, the real values A„ are eigenvalues of T. Here we propose a way 
to construct such operator by using SUSY-QM arguments. 

One of us J2|, was able to consider a p-adic stochastic process having an 
underlying hidden Parisi-Sourlas supersymmetry, as the effective motion of a 
particle in a potential which can be expanded in terms of an infinite collection 
of p-adic harmonic oscillators with fundamental (Wick-rotated imaginary) fre- 
quencies Lu p = i In p (p is a prime) and whose harmonics are u) Pjn = * m P n - 
This p-adic harmonic oscillator potential allowed to determine a one-to-one cor- 
respondence between the amplitudes of oscillations a n (and phases) with the 
imaginary parts of the zeroes of zeta, A„, after solving a inverse eigenvalue 
problem. 

In SUSY-QM two iso-spectral operators H^ and H^ are defined in terms 
of the so called SUSY-QM potential. Here we use the SUSY-QM model proposed 
in J! based on the pioneering work of B. Julia ]2q] , where the zeta-function and 
its fermionic version were related to the partition function of a system of p-adic 
oscillators in thermal equilibrium at a temperature T. The fermionic zeta- 
function has zeroes at the same positions of the ordinary Riemann function plus 
a zero at 1/2 + «0, this zero is associated to the SUSY ground state. See also 



the reference 26 



We propose an ansatz for the following antisymmctrized SUSY QM potential 
(a "p-adic Fourier expansion"): 

#(*; a, b) = n^(«|V + bfp-n nE(4V jI + bfpn, « 

p j p j 



p = primes, j = naturals, a = {a- } and b = {b p }. This comes from the 
following SUSY Schrodinger equation associated with the H^ Hamiltonian 

{jk + $ ) {~L + $ ) ^ +)(x) = A « +) ^ +)(x) ' (2) 

where we set Ti = Ira = 1. SUSY imposes that <&(a;; a, 6) is antisymmetric in x 
so it must vanish at the origin. Hence, <& 2 (a;;a, b) is an even function of x so 
the left/right turning points obey: x£ = —X^ for all orbits, for each n = 1, 
2, ... We define x n = x^'. 

The quantization conditions using the fermionic phase path integral calcu- 
lation (which is not the same as the WKB and sometimes is called the Comtet, 
Bandrauk and Campbell formula pi) are, set h = 2m = 1, so all quantities are 
written in dimensionless variables for simplicity, 

I n (x n ;K;a,b) = 4/ dx [A„ - $ 2 (x;a, b)] = irn, (3) 

Jo 

for n — 1, 2, ... and A n are imaginary parts of the zeroes of zeta. 

The second set of equations are given by the definition of the turning points 
of the bound state orbits: 

$ 2 (x„;a,6) =A„; n= 1,2,... (4) 

So, from the two sets of equations (g) and (|j) we get what we are looking 
for: Amplitudes of p-adic harmonic oscillators, (a, b), and (right) turning points 
x n , depending on all the A„. 

It is very plausible that due to some hidden symmetry of the inverse scat- 
tering problem there may be many solutions for the amplitudes of the p-adic 
harmonic oscillators and turning points; i.e. many different SUSY potentials 
$(x; a, b) do the job for many different sets (a, b) and x n . Numbers related by 
a symmetry which leaves fixed the eigenvalues of the SUSY QM model = imag- 
inary parts of the zeroes of zeta (fixed points). We are not concerned with this 
case now only with proving that a solution exists. 

We emphasize that the integral quantization condition for the energy of the 
orbitals of the SUSY QM model based on the fermionic phase space approxi- 
mation to the path integral Eq. (p|) are only valid for shape-invariant potentials 
V±(x) of the partner Hamiltonians H±. Only for those cases one has that 
the energy levels obtained from this approximation happen to be exact. Since 
imposing shape-invariance amount to an additional constraint one can just for- 
get about this restriction and write down the integral conditions based on a 
fermionic path integral calculation that yields nice approximate results. It is 
far easier to solve the inverse eigenvalue problem in this way than to try to in- 
vert the Shrodinger equations! Therefore, since we are working in the opposite 
direction, and we do not wish to impose an impossible constraint on the SUSY 
potential, and we impose that the energy levels are the imaginary parts of the 



zeroes, then one concludes that the amplitudes of the p-adic harmonic oscilla- 
tor SUSY QM potential, Eq. (pi), (a, b) wee only approximate by our method. 
The numerical results will be given in a forthcoming publication. In the next 
section we propose an approach for a direct proof of the RH. It is based on 
the idea of relating the non-trivial zeroes of the C with orthogonalities between 
eigenfunctions of a conveniently chosen operator. See Q, Q and J22|. 

3 A way to prove the Riemann conjecture 

The essence of our proposal is based in finding the appropriate D\ operator 

D _1_ _^K_ k 

dint dint 

whose eigenvalues s are complex-valued, and its eigenfunctions are given by 

^ s (t)=t~ s+k e v( - t \ (6) 

Notice that D\ is not self adjoint with eigenvalues given by complex valued 
numbers s. 

Also we define a partner operator of D\ as follows, 

d dv 

D > = dnTt + dnTt +k - (7) 

For this operator we have 

D 2 ^ s (t) = (-s + 2k + 2t^)^ s (t). (8) 

The key of our approach relies in choosing the V to be related to the Bernoulli 
string spectral counting function, given by a Jacobi theta series, 



e 2V{t) = 



£ e -™ v 0) 



The Jacobi's theta series is deeply connected to the statistics of Brownian motion 
and integral representations of the Riemann zeta- function J23j . 

Then we are considering a family of D\ operators, each characterized by two 
real numbers k and I to be chosen conveniently. Notice also that D\ is invariant 
under scale transformations of t and F = e v , due to dV/(dlni) = din F/ (dint) 
[ fl6| . The D\ defined in |L6| has k — and a different definition of F. 

Let's recall the functional equation of the Riemann zeta-function |15| , 

Z(s) = n-^T (|) COO - „-W*r (i^) C(l - a) = Z(\ a). (10) 
We define the inner product as follows: 

oo 

(f\g) = jr9j. (ii) 

o 



With this definition, the inner product of two D\ eigenfunctions is 



(V'silV's 



= a I e 2V t~ sl2+2k - 1 dt= ^Z 



a 

1 



-{2k-s 12 ) 



(12) 



where we have defined s±2 = s* + S2 = X\ + 2; 2 + i(j/2 — 2/1 ) and used the 
expressions (||) and (10). a is a constant to be conveniently chosen so that 
the inner product in the critical domain is semi positive definite. The measure 
of integration dh\t is also scale invariant. The integral is performed after the 



change of variables t 



which gives dt/t = (l/l)dx/x, and using the result 



of equation (|13|), given by Voronin and Karatsuba's book |lq| . 

We recall that Z is the fundamental Riemann function, expressed in terms 
of the Jacobi theta series, ui(x) — X)i° exp(— 7m 2 a;) (see Karatsuba and Voronin 
PI), 



n=l 



x x s/2 - l dx 



x s/2 ~ 1 uj(x)dx 



sis - 1 



r s/2-l +x (l-s)/2-l- 



Lo(x)dx 



= Z(s) = Z{l-s). 

(13) 
The right side is defined for all s, i.e., this formula gives the analytic continuation 
of the function Z(s) onto the entire complex s-plane |19f| . 

Having the relation (a/l)Z[(2/l)(2k - s)] = (a/l)Z[l - 2/l(2k - s)], from 
Z[(2/l)(2k — s)] we can obtain the other expression simply by taking k — > k—l/4, 
I — » — I and a —> —a. If, and only if,8fc — 4 = I, then the above discrete 
transformations become: k — > 1 — k, I — > — Z and a — * —a. The importance of 
the relation 8k — 4 = / will be seen shortly. 

From (121) we obtain the "norm" of any state characterized by s = x + iy, 



a 



{i> s \i>s) = '-jZ 



(k-x) 



(14) 



The "norm" of the vectors is the same for all s having the same x. We will 
choose the domain of the values of s, 1 — s such that they fall inside the critical 
domain: < IZe(s) < 1 and < 7Ze(l — s) < 1. We will see soon why this is 
crucial. Recall that at the boundaries, s = and s ~ 1 we have that Z(s) = oo. 

We exclude IZe(s) — and IZe(s) — 1 because Vallee de la Poussin-Hadamard 
theorem says there are no zeroes of £ at x = and x = 1 py] . 

In particular, for the critical line x = 1/2, the value of the "norm" is 
(a/l)Z[2(2k— l)/l]. Since we will choose that I = 4(2fc — 1), the "norm" becomes 



(a/l)Z(l/2) — (a/l)(— 3.97...), for all states in the critical line, independently 
of the chosen value of k. This forces the value of a to be negative, in the critical 
domain, and for this reason we shall fix it to be equal to —I. 

A recently published report by Elizalde, Moretti and Zerbini pl| , which 
contains some comments about the first version of this paper, considers in detail 
the consequences of equation ( O ) . One of them is that equation (|l^) lose its 
original meaning as a scalar product. Despite of this, we will loosely referring 
this map as a scalar product. The Hilbert space inner product property is not 
required so that the eigenvalues can be also negative. The states have real norm 
squared, which need not however be positive definite. Hermiticity requirement 
implies that the states are orthogonal to the reference state and correspond to 
the zeros at the critical line. The problem is whether there could be also zeros 
outside the critical line but inside the critical strip. 

Also we must caution the reader that our arguments do not rely on the valid- 



ity of the zeta- function regularization procedure 1 20 nor the analytic extension 
of the C, which precludes a rigorous interpretation of the right hand side of (O) 
like a scalar product. We simply can replace the expression "scalar product of 
ip si and -0s 2 " by the map S defined as 

S: C®C -» C 

Oi,s 2 ) h-> b{s 1 ,s 2 ) =—Z 

In other words, our arguments do not rely on an evaluation of the integral 
{ip Sl \ip S2 }, but only in the mapping S(si,S2). 

Now we describe our proposal to prove the RH. From now, we shall set k = 1, 
1 = 4. 

Th. 1. If s = x + iX n , where < x < 1 and A n is such that £(l/2 + iA n ) = 0, 
then the states tp s and ipi- s are orthogonal. 

Proof: From (|l||) it follows that 

(^i_.) = -Z(l/2 + tAn) = 0. (16) 

Figure [j] represents those orthogonal states. 

Th. 2. Any pair of states si and s 2 symmetrically localized with respect to 
the vertical line x = 1/2, and such that y\ — yi = 2A n , are orthogonal. 

Proof: We can always find an s = x+iy such that ^ Sl \ip S2 ) = (V'sl^/'i-s) = 0. 
This follows straightforwardly, by using (|12|), sl + S2 = s* + l — s = l~2iy. Then, 
xi + ^2 = 1 and 2/2 — yi = ~ 2y. By equating the arguments of the Z, it follows, 
from theorem 1, that, (ip s \tpi- s ) = only if s — x + i\ n . Therefore, y = A„, 
x\ + X2 = 1 and y\ — yi = 2A n , as we intended to prove. Figure g represents 
those states and we see that for any point x + iy of the complex plane within the 
critical strip, there is a doubly infinity family of orthogonal states to it, given 
by 1 — x + i(y ± 2A TO ) where C(l/2 + A m ) =0. As a special case we have that 
the orthogonal states to the reference state 1/2 + iO are located in the critical 
line. 



Th. 3. The scalar product of two arbitrary states s\ and s 2 within the 
critical strip is the same as the scalar product of a third state S3 of the critical 
strip with the reference state 1/2 + iO. 

Proof: The statement follows directly from the definition (O) written for 
the pairs Si and s 2 , and 1/2+iO and s 3 . By using (ip ai \%j) S2 ) = (^o+iolVVf+aa) = 
('0i/2+io|V ; sJ+s 2 -i/2> I yields directly s 3 = sj + s 2 - 1/2. 

Let's recall that Z(s) = 7r- s / 2 r(s/2)C(s) and (V> Sl |V> S2 ) = ~ z ( l - s i/ 2 ~ 
S2/2) for the chosen values of k and I. This scalar product can be rewritten as 
(ipi/2+io\^s'+s 2 -i/2)- If we define s = 1 - s*/2 - s 2 /2, so that sj + s 2 - 1/2 = 
3/2 — 2s, then the scalar product between any states in the critical domain can 
be rewritten as, 

(ipsi \lpe2) = <^l/2+io|^3/2-2s) = -^( s )- (17) 

Due to the fact that Z(s) = if and only if £(s) = 0, therefore £(s) = if 
and only if "03/2-2s is orthogonal to "01/2+io- Then the RH is equivalent to the 
following statement: The orthogonal states to the reference state are l/2±2iA n . 
Now let's see an interesting consequence of theorem 1. If we define the 
superposition 

|tf> = |V.) + |Vi_.}, (18) 

and evaluate the "norm" of ty, we note that the interference terms are real 
valued 

(ips\i>i- s ) + (1>i-.\4>.) = -Z Q + iy) - Z Q - iyj . (19) 

The arguments of the Z's are in the critical line. If y — X n , for all these states, 
the norm-squared of the sum equals the sum of the norm-squares. This is more 
precise, like the Pythagoras rule. See 116). Therefore, the destructive interfering 
states are distributed on infinite horizontal lines, each of them labelled by one 
of the zeroes of the (. The interference is destructive when y in s = 1/2 + iy is 
such that y — ±A„, with A„ one of the zeroes of the £. 

Since on the critical line the norm-squared given by ( |14| ) is the same for all 
the states labeled by 1/2 ± iy, we can visualize all those states as living on a 
sphere and having the minimal "norm" . 

Notice the form of (|l|). We get as we should a quantity plus its complex 
conjugate for the interference which gives us a real valued number. This makes 
perfect sense because the second inner product must be the complex conjugate 
of the first inner product since we have reversed the order of the inner product. 

Summing up, the interference of the tp s with V'l-s is destructive if and only 
if the values of y n = X n . In such case the norm-squared of the superposition 
equals the sum of the norm-squares of its constituents, this is roughly speak- 
ing Pythagoras theorem. For Pythagorean rational phases and the Riemann 
conjecture (see ]17[|). 

Now we explore some consequences of the mapping given by equation ( |l5| ) . 
Let's suppose that s = x + iy is a generic state orthogonal to the reference state. 
Equation ([l2]) yields, 

(■>Pi/2 + io\ipa) - -Z(3/4 - s/2) = -Z(l/4+ s/2), (20) 



together with their complex conjugates. It follows that if x + iy is orthogonal 
to 1/2 + iO, also orthogonal to this state are: 1 — x + iy, 1 — x — iy and x — iy. 

The arguments of the Z in each single one of those orthogonality relations 
must be of the form 3/4 — x/2 ± iy/2 and 1/4 + x/2 ± iy/2. From the definition 
of Z, there exists a zero if and only if £(x' + iy') = 0. We are going to determine 
relations between s' = x' + iy' and s = x + iy such that the symmetries of the 
Z and the orthogonality condition of states to the reference state are satisfied. 

From the properties of the Ricmann ^-function it follows that if s' is a zero, 
then there are three another zeroes, located at a'*, 1 — s' and 1 — s'*. See figure 
|3|. The 4 arguments of the Z appearing in the orthogonality conditions must 
be related to the location of the 4 zeroes of the £. There are 24 orthogonality 
pairwise combinations, easily identified if we maintain fixed the inner rectangle 
of figure g and perform all permutations of the labels of the orthogonal states. 

If one studies all the possible families of mappings due to all the possible 
4! = 24 permutations of the vertices of the outer rectangle one will inevitably 
introduce constraints among the four vertices of the same rectangle; i. e. like 
s* = 1 — s, and this will trivially lead to the RH. This is unacceptable. 

There are 4 relations of s with s' of the form s = —1/2 + 2x' ± 2iy' or 
s = 3/2 — 2x' ± 2iy', that we call generic cases. They are both reflection- 
symmetric of each other. These simply state a correspondence between the 
position of the orthogonal states and hypothetical zeroes located in any point 
of the critical strip. 

4 of the relations lead to an identification of the vertices of the rectangle 
symmetrically located respect to the critical line, that is s' = 1 — s'* or s'* — 
1 — s', which implies the trivial result that those vertices are located on the 
critical line and the imaginary parts are related by y — ±2y'. These 8 cases are 
compatible with the existence of zeroes in the critical line. 

The remaining 12 relations lead to an identification of all the four vertices of 
the rectangle at the point s' = 1/2+iO which correspond to the state s = 1/2+iO. 
These cases do not correspond to any orthogonality relation, due to in fact 
represent the inner product of the reference state with itself. 

It can not be discarded the presence of orthogonal states to the reference 
located outside the critical line, the generic cases, and in consequence the RH 
can not be deduced from this analysis. 

Notice that the "norm" of the state corresponding to s = is infinite and 
that the states orthogonal to the s = are those states whose s = 1 + iy n — 
1 + 2i\ n . The steps to this proof of the RH did not rely on the same reductio ad 
absurdum argument proposed by Pitkanen (|16)). In our approach, the RH is a 
consequence of the symmetries of the orthogonal states as we intend to prove 
next, then Th. 4 implies a proof of the RH assuming these symmetries. 

Th. 4 ■ The symmetries of the orthogonal states shown in figure 3 are 
preserved for any map S, equation ([l5|), which give rise to Z(as + b), if a and b 
are such that 2b + a — 1. 

Proof: If s = x + iy is orthogonal to a reference state, then the Riemann 
zeta has zeroes at s' = x' + iy' , s'*, 1 — s' and 1 — s'*. If we equate as + b = s', 
then also as* + b = s'*. Now, 1 — s' can be equated to a(l — s) + b and 1 — a'* 



can be equated to a(l — s*) + b if and only if if 2b + a = 1. 

Notice that our election a = —2/1, b = (4k — l)/l is compatible with this 
symmetry if k and I are related by I = 4(2fc — 1). Reciprocally, if we assume 
that the orthogonal states have the symmetries of figure 3, then a and b must 
be related by 26 + a = 1, which gives rise to a very specific relation between k 
and I, obtained from a + 2b — 1; a, b real. The two generic cases correspond 
to taking k = 1, I = 4 and k = 0, I = —4 respectively. It is clear that a map 
with arbitrary values of a and b docs not preserve the mentioned symmetries. 
Theorem 4 contains a genuine proof of the RH assuming that the invariancc 
under the double symmetry is true. 

We have a family of D\ operators, each labelled by (k, I). Since we can set 

I = 4(2fc — 1) due to the constraint 1 = a(k, I) + 2b(k, I) imposed by the double- 
ts) 
reflection symmetry, we can parametrize the eigenfunctions as ip. s , where k 

and s are continuous variables. Let's consider two of those operators, D[ 1 ' and 
D{ . It must exist a one to one correspondence between eigenfunctions of this 
pair of operators at any given point s. We can see that to prove the RH is 
equivalent to the following statement: 

St. (i). If t = 1/2 + «0 is the reference states and if si and s 2 are both 
points having x\ = a; 2 = x, then if V'si is orthogonal to %j)\ , then ips 2 2 is 
also orthogonal to i)r t . In other words, the orthogonality of states (with the 
same x) to the reference state is independent of / and k. 

Due to the fact that the inner products of the states of (i) with the corre- 
sponding reference state are given by Z[(2/li)(2ki —t* — S\)} and Z[(2/l2)(2k2 — 
t* — S2)], then the orthogonality gives rise to two zeroes of the Riemann zeta. 
Those zeroes are of the form s[ = x' + iy[ and s' 2 = x' + iy 2 , that is, they are 
to be located on the same vertical line. 



This immediately allows us to write s[ — (2/li)(2ki — t* — si) and 



2 



(2// 2 )(2fc 2 — t* — S2), whose real parts are respectively x[ = (2/li)(2ki — 1/2 — x) 
and x 2 = (2// 2 )(2fc 2 - 1/2 - x). If, and only if, x = 1/2 and h = 4(2fci - 1), 
Z 2 = 4(2fc 2 - 1), then it follows that x[ = x' 2 = x = 1/2, y[ = (—2/h)yi, 
y' 2 = (— 2/? 2 )y 2 . Notice that the values of j/i, y 2 do depend on (k,l). y[ = \\ 
and y' 2 = A 2 are the imaginary parts of the two zeroes of £. Then (i) is equivalent 
to the RH. 

These arguments can be generalized for all (k, I) obeying I = 4(2fc — 1) 
as follows. From s' = a(k, l)s + b(k, I) and a + 2b = 1 we can write s' = 
a(s — 1/2) + 1/2. So x' = 1/2 + a(x — 1/2) and, for a fixed x the position of 
a zero could continuously change by changing (fc, I). If the zeroes are supposed 
to form a discrete set, then one must have that x = 1/2, so that x' = 1/2 is the 
only consistent value one can have. 

Another way of rephrasing this is saying that the family of the D[ ' op- 
erators yields a continuous family of pairs (x,x') having the double-reflection 
symmetry x — > 1 — x, x' — > 1 — x', from which one arrives at the constraint 1/2 = 
a(k, l)/2 + b(k, I). This means that all the curves given by x' = a(k, l)x + b(k, I 
must contain the common point in the x — x' plane given by (1/2, 1/2), for 



all values of (k,l), obeying I = 4(2fc — 1). Also, if ips is orthogonal to the 
reference state 1/2 + iO then s' is a zero of the £, and the real parts are related 
by x' = a(x — 1/2) + 1/2. Due to the statement (i) this real part must to be 
independent of k namely, independent of a. This can be satisfied only if the 
orthogonal state satisfies x = 1/2, which gives x' — 1/2, the RH. 

Let us summarize what we have found so far, for < IZe(s) < 1. When we 
choose k = 1, I = 4 then we have the following results (i), (ii), (iii). (i) The 
interference between ip s and V'l-s is destructive at the horizontal lines defined 
by < x < 1 and y = ±A„. (ii) If two states Si = 1/2 + iy\ and s 2 = 1/2 + iy 2 
on the critical line are orthogonal, then y\—yi = ±2A ra for one given n. (iii) Any 
pair of orthogonal states must obey X\ + X% = 1 and y% — y\ = 2A„. And when 
we assume the double reflection symmetry, then for all (k, I) obeying 8k — 4 = I 
we have: (iv) The RH is a consequence of the assumption that the orthogonal 
states to the reference state have the same symmetry like the zeroes of £ have. 

4 The Riemann fractal string and its dual 

Finally, we recall that the e 2V defined by (ph, expressed in terms of a; = i 4 , 
is nothing else than the energy partition function of a Bernoulli string (whose 
standing waves have for wave vectors integral multiples of the inverse of length 

II)- 

Several new relations of Riemann zeta-function with the spectral proper- 
ties of different physical systems have been found during the last years. The 
authors discuss the Epstein zeta-functions [G~3| when dealing with the Laplace 
operator in square/rectangular domains. The counting function of the spectral 
eigenvalues of 2D-Laplacian gives the Epstein zeta. Those functions also ap- 
pear in the physics of p-branes moving in hyperbolic spaces (negative constant 
curvature); i.e. in the calculation of the effective potential of toroidal p-branes 
living/embedded in target hyperbolic spaces [fl4| 

The spectrum of toroidal p-branes and branes moving in hyperbolic spaces 
of constant negative curvature, is linked to the Epstein zeta-functions [ [l3[ , 

Ce^( S )=E K+nl+ ^ +i)s/2J (21) 

where {n} = {ni,ri2, ...n p +i} is a non-zero integer vector. 

The Riemann zeta is a special case of the Epstein zeta and has a relation 
to the spectrum of fractal strings (Alain Connes had urged to use a different 
word to avoid confusion with string theory, suggests call them "fractal harps" 
instead). Two functions describe the fractal string, the geometric length count- 
ing function Zl{s) and the frequency counting function Zp(s) that are related 
by the fundamental relation involving the Riemann zeta, 

Z F {s) = Z L (s)C(s). (22) 

Let's imagine a chain made out of many links of sizes 1%, I2, ■■■ lj, ■■■ The 
natural frequencies associated with the oscillations of each one of those links 
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are given by: l/h; l/h, ■■■ !/(;> ■■■ and the excitation states are just integer 
multiples of 1/lj'. k/lj] for k = 1, 2, ... j, ... There are described by the geometric 
length counting function Zl(s) and the frequency counting function Zp(s). The 
boundaries of the fractal strings are Cantor sets. See [E] and R. 

Since we have selected the V(t) of Eq. (^) to be related to the Bernoulli string 
spectral counting function we must continue with the fractal strings duality 
discussion. It is not a coincidence that choosing V(i) to be related to the 
Bernoulli string was telling us something! 

The length counting function is the sum, 



oo 



z l (s)=J2 1 j- ( 23 ) 

J'=l 

The frequency counting function is the double sum, 

oo oo 

fc=U=l 

And we always have Zp(s) = Cgen(s) Z 'l(s) . In the case of a fractal string it 
reduces to Eq. (p2|) . The book is entirely based in choosing particular examples 
of fractals which fixed the lengths lj a priori which enable to evaluate the sums 
explicitly. The authors of p3| were able to find the complex dimensions of the 
fractal strings by finding the location of the poles of the Zl(s). It was essential 
that the poles do not coincide with the zeroes of C( s )- 

Let us explain how one obtains the complex dimensions of the Cantor string. 
The sequence of lengths must always be decreasing: l± > 1% > ... or at most 
equal. Lets take the Cantor set as an example; there is 1 segment of length 
equal to 1/3; 2 segments of length = 1/9; ... 2™ segments of l n = (1/3)™ and so 
forth. The degeneracy of each link of length l n is w n — 2™. So the geometric 
counting function of the Cantor string is: 

/1\™ S 1 

^(*) = £(<«)* = £ 2n ( 3) =TT2^- (25) 

n n ^ ' 

so the poles of this sum yield the complex dimension: 1 = e l2rlTr = 2.(3) _s ; 
taking the logarithm on both sides yields s„ — (ln2/lnZ) ± i(2n7r/ln3) for the 
complex dimensions for n — 0, 1, 2, ... 

The real part is the standard dimension of the Cantor set. The imaginary 
parts are periodic whose period is 27r/ln3. 

We propose that the physical meaning of these complex valued dimensions 
may be relevant for a theory of quantum gravity, because of interference of 
complex dimensions. 

It is very important to insist on the condition < 5Rs < 1 since a fractal 
string of real valued dimension D, embedded in a space of R d , must satisfy the 
constraint: d — 1 < D < d. 
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The Cantor string, the golden string and the Fibonacci string are defined 
respectively for the finite number of scaling ratios given by 

n=r 2 = 1/3; n = 1/2, r 2 = l/2 1+< ^, <f> = 0.618...; n = 1/2, r 2 = 1/3. (26) 

We shall construct the Riemann fractal string (RFS) by different procedures 
than the authors |l3| , as the fractal boundaries of the open 2-D domain. 

The duality relation of the Riemann zeta Z(s) — Z(\ — s) only makes sense 
if, and only if, we write the continuum limit for the length counting function 

/>oo />oo 

Z L {s)= l(x) s dfi(x) andZ L (l-s) = / f(a:) 1 - s d/I(a:), (27) 

Jo Jo 

and the frequency counting function 

/>oo />oo 

Zf{s) = / f{x)~ s dn(x) and Z F {\ - s) = / f(x)- 1+s dfi(x), (28) 

Jo Jo 

where (l(x), l(x)), {f{x), f(x)), (/J>(x), fi(x)) are the respective complex- valued, 
and their complex conjugates, lengths, frequencies and measures associated with 
the continuum limit of the RFS and their dual, that in the general case are 
complex valued maps from the real line to the complex plane. The measure 
fi(x) is such that J / °° d/j(x)\ = 1. 

In reference |g| a detailed analysis of the vector calculus and contour inte- 
grals on fractal curves (boundaries of a bounded domain in the complex plane) 
and interfaces is given. This was attained by constructing pseudo-measures of 
integration based on iterated function systems. Expressions for the contour 
integrals are obtained by means of a suitable renormalization procedure (the 
length of a fractal contour is infinite) and the solution of the Dirichlet problem 
on bounded two-dimensional domains possessing fractal boundaries is given. 
These analytical tools allow us to find in principle the complex valued functions 
appearing in equation ( |29| ) and to define the integrals in ( p7||2^ ). 

From (g7p and ( p8[) one can obtain the generalization of the Lapidus and 
Frankenhuysen result, our proposal for RFS, 



z F = f^ttxr*Mx) 

Zl J °° l(x) s d/j,(x) 



CW, (29) 



with analogous relation for the dual string. Lapidus and Frankenhuysen results 
are recovered by using a Dirac delta distribution for the measure. 

When s is one of the zeroes of the £, the equation (E9h implies that the 
geometric length counting function has a simple pole or that the frequency 
counting function has a simple zero. On basis of this, we define the Riemann 
fractal string as such string whose geometric length counting function has a 
simple pole precisely at the zeroes of the Riemann zeta-function. 
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5 Concluding remarks 

In a previous work [EJ we found very suggestive relations of the golden mean and 
the distribution of the imaginary parts of the zeroes of zeta. We pointed out 
that this could be related to the multifractal character of the prime numbers 
distribution (See |p5| ). Here we observed yet another evidence of the multi- 
fractal distribution of primes in the similarities in Wolf's formulae with those 
of Lapidus and Frankcnhuysen: 

p{Si)~U, p(Si) q ~lt, q~s, Xq {l) ~ Z L {s), l~L, (30) 

where L is the initial length where one begins to construct the fractal string by 
defining N finite segments of lengths lj — rjL, j — 1,2, ...N so that what is 
left is L(l — R) where R = £\ rj < 1. One chooses L so that L(l — R) = 1 
which means that the first length of the iteration has unit length. From Wolf's 
paper we have that p(Si) is the measure of the set Si, the Xq are the moments 
or partition functions, which scale like a l T ( q \ 

The Riemann zeta is deeply connected to fractal strings (whose boundaries 
are Cantor sets) and it is no wonder why it is related to quantum chaos, random 
matrix models, random walks, Brownian motion, etc. 
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Figure 1: Lines between dots and crosses represent pairwise orthogonal states. 
Dots represent zeroes of the £. Crosses on the critical line represent states 
orthogonal to the 1/2+iO. Crosses on the IZe(s) = 1 represent states orthogonal 
to the + i0. Crosses on the IZe(s) = represent states orthogonal to the 1 + iO. 
The states 1/2 + i\ n arc orthogonal to 1/2 — i\ n . On the critical line, pairs of 
dots and/or pairs of crosses are mutually orthogonal. Notice that for simplicity 
we are representing the orthogonalities of states having only x = 0, x = 1/2 and 
x = 1. Here we are referring the particular case k = 1, I = 4. 
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Figure 2: The state associated to a point x + iy is orthogonal to a doubly 
infinite series of states at 1 — x + i(y ± 2A„), for n = 1,2, ...oo. Here we are 
referring the particular case k = 1, 1 = 4. 
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Figure 3: The dots represent generic zeroes of the £. The crosses repre- 
sent generic states orthogonal to the reference state 1/2 + iO. The numbers 
3/4 — x/2 — iy/2, etc, are the arguments of Z appearing in the orthogonality 
relations between states orthogonal to the reference state. Due to the functional 
equation of the Riemann zeta- function (J10|) , these arguments are just the aver- 
age values between 1/2 + iO and those orthogonal states. Here we are referring 
the particular case k = 1, 1 = 4. 
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